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NEAR FIELD TRAILING EDGE TONE NOISE COMPUTATION 
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Taitech, Inc. 

Cleveland, Ohio 44135, USA 
e-mail: fsloh@turbot.grc.nasa.gov 


Abstract 

Blunt trailing edges in a flow often gener- 
ate tone noise due to wall-jet shear layer and 
vortex shedding. In this paper, the space- 
time conservation element and solution ele- 
ment (CE/SE) method is employed to numeri- 
cally study the near-field noise of blunt trailing 
edges. Two typical cases, namely, flow past a 
circular cylinder (aeolian noise problem ) and 
flow past a flat plate of finite thickness are con- 
sidered. The computed frequencies and the es- 
timated sound pressure level (SPL) compare 
well with experimental data. For the aeolian 
noise problem, comparisons with the results of 
other numerical approaches are also presented. 

1 Introduction 

In today’s aviation industry, everyday, thousands of 
commercial airliners take off and land at the airports 
around the world, accompanied by various uncomfort- 
able noises, polluting the environment of the nearby ar- 
eas. The noises due to airframe components, e.g. airfoil 
slat, landing gear and opeings, etc. are one of the major 
noise sources. They are generated by the interaction of 
the vortex street in the turbulent wake. A blunt trailing 
edge often generates certain tone noise in the flow. Nu- 
merical investigations of airframe noise based on direct 
numerical simulation (DNS) need very fine grids and are 
hence expensive. Large eddy simulations (LES) allow 
a coarser grid since only large eddies are resolved and 
the smaller ones are modeled by a subgrid scale (SGS) 
model. Furthermore, some researchers use Reynolds 
averaged Navier- Stokes (RANS) methods to save more 
storage and CPU time [1-5]. 

In the present paper, we use a DNS type approach to 
investigate two benchmark noise problems, namely, the 
noise generated by the flow past a circular cylinder (ae- 
olian noise) and the noise generated by flow past a flat 
plate of finite thickness. Both may be categorized as 
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blunt trailing edge noise. By ‘DNS type’, we mean that 
the grid may not be fine enough for the Kolmogorov scale 
of turbulence but is fine enough for the large turbulence 
structure (or large eddies) associated with the tone noise 
frequencies. This approach seems somewhat similar to 
LES but there is no explicit filtering or SGS model. More 
details are given later. 

The space-time conservation element and solution el- 
ement method (CE/SE) [6-8] is adopted for the com- 
putation. As demonstrated in our previous papers, the 
CE/SE scheme is well suited for aeroacoustics computa- 
tion [9,10]. Because of the CE/SE non-reflecting bound- 
ary conditions (NRBC), which are based on flux balance 
[11], a smaller near field computational domain can be 
used in the present numerical simulation and helps to 
save both memory and CPU time. 

The 2-D unstructured-grid Navier- Stokes (N-S) 
CE/SE scheme used here is briefly discussed in Section 
2. Section 3 illustrates the initial and boundary condi- 
tions for these blunt trailing edge noise problems, in par- 
ticular, the CE/SE NRBC. The numerical results are 
presented and compared to available experimental 
findings [12,15] in Section 4. Conclusions are drawn 
in Section 5. 

2 The Numerical Scheme 

The space-time conservation element and solution ele- 
ment (CE/SE) method is a low order (second order accu- 
rate in time and space) scheme with low dissipation. It is 
chosen in the present numerical investigation. A detailed 
description of the method can be found in [6- 8] . In this 
Section, only a brief sketch of the scheme is given. 

2.1 Conservation Form of the 2-D Unsteady 
Compressible Navier-Stokes Equations 

In general, the CE/SE method systematically solves a 
set of integral equations derived directly from the physi- 
cal conservation laws and naturally captures shocks and 
other discontinuities in the flow. Both conservative vari- 
ables and their derivatives are computed simultaneously 
as unknowns. 
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Consider a dimensionless conservation form of the un- 
steady 2-D Navier- Stokes equations of a perfect gas. Let 
p, u , v, p , and 7 be the density, streamwise velocity com- 
ponent, transverse velocity component, static pressure, 
and constant specific heat ratio, respectively. The 2-D 
Navier-Stokes equations then can be written in the fol- 
lowing vector form: 


Ut + F x + G y — 0, (1) 


where x, y, and t are the streamwise and transverse co- 
ordinates and time, respectively. The conservative flow 
variable vector U and the flux vectors in the streamwise 
and radial directions, F and G, are given by: 


U = 


{Ui 


u 2 

Us ’ 

\u 4 J 


/ F1 \ 

F 2 

F 3 

\F 4 J 



with 

Ui = p, U 2 = pu , Us = pv, 

U 4 =p /( 7 - 1) + p(u 2 + v 2 )/2. 

The flux vectors are further split into inviscid and viscous 
fluxes: 

F = Fi — F v ,G = Gi-G v , 

where the inviscid fluxes are the same as in the Euler 
equations: 

F n = U 2j 

F i2 = (7 - 1)^4 + [(3 - 7 Wi - (7 - 1 Wl] /2 U u 
F i3 = U 2 U 3 /U U 

Fi 4 = 7^2t/ 4 /t/i - (7 - 1 ) 1^2 [Ui + t/ 3 2 ] / 21 /f, 

Gii = U 3 , Gi 2 = I/2E/3M, 

Gi 3 = (7 - 1 )U A + [(3 - 7 )Ul - (7 - 1 )G|] /2I7 X , 
Gi 4 = lUsU A lU x - (7 - l)I/ 3 [Ui + G 2 ] /2f/ 2 , 



hexagon cylinder ADBECFA - A’D’B’E’C’FA’ 
wilhils3CEs:OADB-aAD’B’ 5 OBEC-O’B’E’C’, 
and OCFA - O’C’FA’ 


Figure 1: CE/SE unstructured grid in space-time E 3 
space. 


7 d U 4 u 2 +v 2 

Prdy^lh 2 ^ 

where u, v, u x , u y , , v y are respectively the a:— and p— 

flow velocity components and their derivatives, which 
can be written in terms of the conservative variables 
U\,U 2 , U 3 and C/4, with Pr (= 0.72) being the Prandtl 
number, p the viscosity, the velocity divergence 

V*"V I/-# + Vy. 

By considering (x,y,t) as coordinates of a three- 
dimensional Euclidean space, E 3 , and using Gauss’ di- 
vergence theorem, it follows that Eq. (1) is equivalent to 
the following integral conservation law: 


and the viscous fluxes are: 

F v 1 = 0, F v 2 = n(2 u x - | V-V), 

Fv 3 — P j {v x + 

F v 4 — p[2'U , u^ + ( u y T" v x ^)v ^ (^*"V") < uT _ 
7 9 ,£/ 4 u 2 + u 2 

2 ^ 

G^l — 0, G^2 — f^{V X + Vjy) , 

G^ 3 = - §V-V), 

G ^4 = p[2w y + (u y + v x )u - |(V-V)u+ 


<b H m dS = 0, m = 1, 2 , 3, 4, (2) 

J S(V) 

where 5(F) denotes the surface around a volume F in 
E /3 and H m = {F m ,G m , U m ) . 

2.2 Unstructured Grid for CE/SE 

The CE/SE scheme is constructed to take advantage of 
an unstructured triangle grid. The unstructured geome- 
try used with the CE/SE scheme is illustrated in Fig. 1. 
Here, AABC is a typical triangle cell and D,E,F are 
the triangle centers of the neighboring cells. The flow 
variables at the previous time step are stored at the trian- 
gle cell centers. Three quadrilateral cylinders (conserva- 
tion elements) are formed by the edges that connect the 


NAS A/CR— 2002-2 12013 


2 

American Institute of Aeronautics and Astronautics 


vertices and the center of the triangle and its three neigh- 
bors. In the space-time E 3 space, Eq. (2) is applied to 
the hexagon cylinder ADBECF — A DBE 1 C'F ' of 
volume V that consists of these 3 quadrilateral cylinder 
CE s, see Fig 1. 

In the CE/SE scheme, the above flux conservation re- 
lation, Eq. (2), in space-time is the only mechanism that 
transfers information between node points. A conserva- 
tion element CE (here, quadrilateral cylinders) is the fi- 
nite volume to which Eq. (2) is applied. Discontinuities 
are allowed to occur in a conservation element. A so- 
lution element SE associated with a grid node (e.g., D, 
E , or F in Fig. 1) is here a set of interface planes in E 3 
that passes through this node (e.g. DAA D 1 , DBB D ', 
EBB E 1 , ECC E etc.). Each surface S(CE) is made 
up of segments belonging to two neighboring CD’s. 
Within a given solution element SE(j,n ), where j, n are 
the node index and time step, respectively, the flow vari- 
ables are not only considered continuous but are also ap- 
proximated by the linear Taylor series expansions. The 
surface flux can be calculated accurately and easily by 
first evaluating the flux vectors at the geometrical cen- 
ter of the surface through the above Taylor series expan- 
sions. 

At time level n, the solution variables U, U x , and 
Uy are given at the three nodes D,E,F in Fig. 1 and 
U, U x and U y at O' at the new time level n + 1 are 
to be computed. In principle, each of the 3 CE s pro- 
vides 4 scalar equations when Eq. (2) is applied to the el- 
ement. Hence, the 12 scalar equations needed for the 12 
scalar unknowns at O f are available. All the unknowns 
are computed based on these relations. No extrapolations 
(interpolations) across a stencil of cells are needed. 

More details about the unstructured CE/SE method 
can be found in [ 8 ]. 

3 The Trailing Edge Noise Problems 

It is well known that when the trailing edge of a body 
is sharp, the turbulent wake generates broadband noise, 
while a blunt trailing edge produces certain tone noises 
due to the presence of wall-jets and the consequent large 
scale turblulence. Typical blunt trailing edges include, 
e.g., the flat plate of finite thickness and the circular 
cylinder. They are investigated separately in this section. 

3.1 The Aeolian Noise Problem 

Aeolian noises, or the sounds generated by flow past a 
circular cylinder, are closely related to airframe noise 
and noise in automobiles , e.g. noise induced by auto- 
antenna. As it is important, the problem is chosen as one 
of the benchmark problems in CAA Workshop II [1]. 

Consider a Mach number M = 0.2 uniform flow past 
a 2-D circular cylinder of diameter D = 1.9 cm. Fig. 2 
shows a sketch of the computational domain and the tri- 
angulated grid. As shown in the figure, the grid is ob- 


tained by dividing a rectangular cell into 4 triangles. By 
so doing, the advantage is that the cell size, which is im- 
portant in aeroacoustics computation, is well-controlled. 
The diameter D of the circular cylinder, the ambient 
speed of sound ao and the ambient flow density are cho- 
sen as the scales of length, velocity and density respec- 
tively. 

top NRBC 



bottom NRBC 

Figure 2: Geometry of the computational domain and 
unstructured grid, —2. ID ^ x ^ 5.5D, —4 D ^ y ^ 
4D; 7 exponentially growing rectangular cells (as sponge 
zones) at the top and bottom are not shown. 


3.1.1 Initial Conditions Initially, the flow of the 
entire domain is set at the ambient conditions, i.e., 

Po = l.o, Wo =0, Vo = 0, Po = 1 / 7 , 

All spatial derivatives are set to zero. 

3.1.2 Non-Reflecting Boundary Condition 
(NRBC) for CE/SE In the CE/SE scheme, NRBCs 
are constructed so as to allow fluxes from the interior 
domain to a boundary CE smoothly exit to the exterior 
of the domain. There are various variants of the NRBCs 
[9-11], the following are the ones employed in this paper. 

For a grid node (j, n) lying at the outer border of the 
domain, where j is the grid node index number and n the 
time step, the Type I NRBC requires that 

(U X )? = ( Uy )? = 0 , 

while ITj is kept fixed at the initially given steady bound- 
ary value. At the outflow boundary, where there are 
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substantial gradients in the radial direction, the Type II 
NRBC requires that 

(■ u x )? = 0 , 

while U'j and (U y )J are now defined by simple extrap- 
olation from the nearest interior node j\ i.e., 

u? = ur 1 ( u v)j = ( u v)t 1 - 

As will be observed later, these NRBCs are robust 
enough to allow a near field computation without dis- 
turbing or distorting the flow and acoustic fields. 

3.1.3 Boundary Conditions At the inflow 
boundary, the flow variables are specified as the M = 0.2 
uniform flow: 

Pi = 1.0, Ui = 0.2, Vi =0, Pi = 1 / 7 , 

and their spatial derivatives are conveniently set to zero. 
Note that these inlet conditions also play a role as ab- 
sorbing boundary conditions in the CE/SE scheme. At 
the surface of the circular cylinder, no- slip wall condition 
is applied. At the domain top and bottom Type I NRBC 
of the CE/SE method is used. The CE/SE outflow NRBC 
(Type II) is applied at the outflow boundary. Due to the 
inconsistency of the boundary conditions at the top, bot- 
tom and the inflow one, unpleasant non-physical gradi- 
ents appear at the upper and lower left corners, although 
they do not propagate to the interior. To remove them, 7 
rectangular cells with exponentially fast growth in size 
in the y direction are added to the top and bottom as 
“sponge zones”. Each of these rectangular cells is then 
divided into 4 triangle cells as in the interior. 

3.2 The Trailing Edge Tone Noise of a Flate 
Plate of Finite Thickness 

In the case of a flat plate of finite thickness, due to flow 
separation, a wall-jet is formed downstream of the blunt 
trailing edge of the flat plate, tone noise is thus generated 
from the consequent Karman vortex street. 

Consider a Mach number M = 0.3 uniform flow past 
a 2-D flat plate of thickness D = 1” (25 4cm). The 
flat plate extends into the computational domain with 
a length of 3D. Fig. 3 shows a sketch of the compu- 
tational domain. Similar to the previous problem, The 
plate thickness D , the ambient speed of sound ao and 
the ambient flow density are used to scale length, veloc- 
ity and density respectively. The domain shown in Fig. 3 
is exactly the computational domain, no sponge zone is 
used for this problem. As in the previous problem, the 
‘unstructured’ grid is obtained based on cutting a struc- 
tured rectangular cell into 4 triangles. 



bottom NRBC 


Figure 3: The computational domain for the flat plate 
noise problem; totally 77,600 triangular cells, no sponge 
zones. 

3.2.1 Initiai Conditions Initially, the flow of the 
entire domain is set at the ambient conditions, i.e., 

Po = l.o, u 0 =0, v 0 = 0, Po = 1 / 7 , 

All spatial derivatives are set to zero. 

3.2.2 Boundary Conditions At the inflow 
boundary, the flow variables are specified as the M = 0.3 
uniform flow: 

pi = 1.0, Ui = 0.3, Vi = 0, Pi = I/ 7 , 

and their spatial derivatives are conveniently set to zero. 
The no- slip wall boundary condition is applied to the flat 
plate. The top, bottom and outflow NRBC are the same 
as the previous problem, but no sponge zones are used. 

4 Numerical Results 

In this section, numerical results for both the aeolian 
and flat plate trailing edge noises are presented and com- 
pared to their respective experimantal data and results by 
other numerical approaches when they are available. 

4.1 Aeolian Noise Problem 

In order to have a sizable sampling time series for FFT 
(fast Fourier transform) analysis, with non-dimensional 
time step size At = .005, over half million time steps 
are run. Fig. 4 - Fig. 6 demonstrate the isobars or numer- 
ical Schlierens at time step 410,000 and 545,000. The 
unsteadiness and the Karman vortex street in the wake 
are clearly displayed. A point located at (0D,3.25D) 
(Fig. 4) is selected for recording the time series. At this 
location, there is no direct inluence from the Karman vor- 
tex street. Fig. 7 illustrates the non-dimensional pressure 
history at this point. Fig. 8 depicts the result of FFT, i.e., 
the sound pressure level (SPF) vs. frequency at this se- 
lected location. The SPF is based on the relative level to 
20 //Pa. 

4.1.1 Tone noise frequency At Reynolds 
number Re= 90, 000 the highest peak of SPF in the 
plot corresponds to a computed frequency of 617 Hz 
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7.7D 


FFT pt. 


8 ) 


|M = 0.2, isobars at time step 410,0001 


15.4 R 


aeolian noise problem (M=0.2) 



'/! 


isobars at time step 545,000 
• minmax: 0.6414, 0.7817 



Figure 6: Instantaneous isobars at time step 545,000. 


Figure 4: Instantaneous isobars at time step 410,000. 



numerical Schlieren 
at time step 410,000 

Figure 5 : Instantaneous numerical Schlierens at time step 
410,000. 


(with a 3 db ‘bandwidth’ or binwidth of 88 Hz), which 
is quite close to the experimental one of 643 Hz [1] at 
this M — 0.2 Mach number, despite the crude resolution 
in frequency. A numerical experiment with a “very 
high” Reynolds number is also tested. The computed 
frequency is 652 Hz, slightly higher but still within a 4% 
range of the experimental 643 Hz. The phenomenon that 
the frequency ( or Strouhal number) around this range of 
Reynolds number does not significantly change with Re 
is justified by experimental observation (see Schlichting 

[13] ). and also agrees with Townsend’s observation 

[14] that the large scale turbulence only weakly depends 
on the Reynolds number. The term “very high” used 
here is a synonym of turning off the viscous terms in 
the scheme. As a result of the numerical damping in 
the scheme, viscous effects still exist. According to the 
chart in Schlichting’s book (p.32, [13] ), we speculate 
that the effective Re is probably in the order of O(10 4 ) 
to O(10 5 ). 

It is interesting to compare the computed frequency 
with the results by other numerical approaches from the 
CAA Workshop II [12]. 

Table 1 lists the computed frequencies by several other 
researchers [2-5]. Generally, they use either compress- 
ible (Comp.) or incompressible (Incomp.) flow solver 
with RANS (Reynolds averaged Navier- Stokes), LES, or 
DNS and then an acoustic solver is applied. It is ob- 
served that, compared to the experimental result of 643 
Hz, there is an error varying from 10 — 60%, while for 
Re from 90, 000 to 0(1 0 5 ), the current method yields a 
error within only 4%, and no acoustic solver is employed 
at all. 
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Table 1 : Comparison of computed frequencies 


Experiment 

Brentner 
et al. 

Kumarasamy 
et al. 

Pope 

Spyropoulos 
and Holmes 

CE/SE 

643 Hz. 

815-1031 

710 

753 

550 

617 

652 

Flow 

Comp. 

Incomp. 

Incomp. 

Incomp. 

Comp. 

Solver 

RANS 

RANS 

DNS 

LES 

DNS 

Acoustic 

Solver 

FW-H 

Helmholtz 

A/VS 

Curie 

none 

Re 

lk-90k 

90k 

20k 

90k 

90k 

very 

high 


pressure history at (-1 .75D, 3.25D) 



Figure 7 : History of pressure fluctuation at the FFT point 
(recorded every 20 time steps). 


SPL at (0.3.25D) in dB 



Figure 8: Sound pressure Level for Re=90,000 at (0, 
3.25D), the highest peak corresponds to 617 Hz and 138 
dB. 
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4.1 .2 Sound pressure level of the tone noise 

Since the computation is two-dimensional, and the com- 
putational domain is in the near field (within 8 diame- 
ters), there is no direct comparison of the sound pressure 
level (SPL) to the experimental data which is of 3-D na- 
ture and measured at a location 35 diameters away from 
the circular cylinder center [12]. However, after some 
approximation and simple analysis from the first princi- 
ple, an estimated SPL for 3-D flow SPL at ro = 35 D in 
the direction of 90° can be obtained. 

Based on the 2-D computation as shown in Fig. 8, at a 
point of (0D, 3.25D) (see the ‘FFT point’ in Fig. 4) which 
is in the 90° direction with x axis, the SPL is 138.5 dB. 
If the noise source is regarded to be a simple one and lo- 
cates approximately at the circular cylinder centerline, 
for 2-D flow, the SPL attenuation is inversely propor- 
tional to r 2 , where r is the distance from the source to 
the measuring point. Therefore, at the point Q (0, 35D) 
in the 90° direction the SPL is estimated to be: 

SPL 2d = 138- 20log(35D/3.25D) = 117.9 dB. 

This result is obtained under the 2-D assumption that the 
circular cylinder is infinitely long. 

According to the experiment [3, 12], a cy liner length 
of 10D is used. Q is assumed at the mid-plane normal 
to the cylinder centerline ( the z direction), the cylinder 
rod ranges from z = —5 D to z — 5 D. For truly 3-D 
flow, let pfy) be the r.m.s. pressure contribution from a 
line segment dz at z (along the cylinder centerline) to Q. 
Then 

r 2 

P® Po=P(0). 

Generally, let L be the half length of the cylinder, the 
square of the total sound pressure at Q is 

£ Po [1 + ( -g / ro ) 2]2 = r 0 pl[a 0 +0.5sin(2a 0 )], 

where ao = tan~ 1 (L/ro). Let L = 5 D and L 
respectively and the ratio 

r = Ezr = 0.1794612 

Poo 

Therefore, the sound pressure level at Q is estimated as: 

SPL 3d = SPL 2d + 10 logR = 110 AdB 

Table 2 lists the numerical SPLs at these two locations at 
an angle of 90° and comparison to the experimental data. 
Recall that the basic computation is a two dimensional 
one, while the experimental data is for a circular cylinder 
of finite length in 3-D space, the agreement is excellent. 


Table 2: SPL of aeolian tone noise at 90° angle 


r 

(distance) 

Experiment 

[12] 

CE/SE 

numerical 

3.25D 


138.5 dB 

35D 

111 dB 

110.4 dB 
(estimated) 


presure history 



time in sec. 


Figure 10: Pressure history at the FFT point (recorded at 
every 2 time steps). 

4.2 Trailing Edge Noise Tone of a flat Plate 

With non-dimensional time step size At = .005, 390,000 
time steps are run. Fig. 9 demonstrates the instantaneous 
isobars at the last time step. The unsymmetric unsteadi- 
ness and the Karman vortex street in the wake are evi- 
dent. Since there is no sponge zone, the ‘matching layer’ 
between the steady Type I NRBC at the top and bottom 
boundaries are observed. Also, at the upper and lower 
left corners, non-physical gradients due to inconsistency 
of inflow boundary condition and the Type I NRBCs are 
present. However, they have practical no influence on the 
interior flow. A point located at (1 3D , 3D) is selected 
for recording the time series. At this location, there is no 
direct inluence from the Karman vortex street. Fig. 10 
illustrates the non-dimensional pressure history at this 
point. Fig. 11 depicts the PSD vs. frequency at this 
selected location and the experimental results by Heine- 
mann et al [15] for similar cases. The Reynolds number 
is in the magnitude of 10 4 to 10 5 (“very high”). The 
highest peak of PSD in the plot corresponds to a com- 
puted frequency of 791 Hz (with a 3 db bin width 33 Hz). 
Experimentally , from Heinemann’s chart in Fig. 11, for 
Mach number M = 0.3, the corresponding Strouhal 
number S = 0.195. Assuming D = 2.54 cm and speed 
of sound equal to 343 m/s, it is equivalent to 790 Hz, 
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FFT data pt. 



Figure 9: Instantaneous isobars at time step 390,000, for flow past the flat plate. 



Experimental 
data for flat 
plate vortex 
street noise: 
790 Hz. 

(Heinemann et 
al.) 


CE/SE result 


791 Hz 



-9.0 -| 1 i i i I I i 

0. 2500. 5000. 7500. 10000. 

Freq. Hz. 


Figure 1 1 : PSD and experimental data. 


Table 3: Flate plate tone noise frequency 


Experiment 
(Heinemann et al) 

CE/SE 

numerical 

790 Hz. 

791 Hz 


almost identical to the computed value, as shown in Ta- 
ble 3. No SPL data is collected since there is no experi- 
mental data in [15] to compare with. 

5 Concluding Remarks 

In this paper, a DNS type numerical approach is 
adopted to simulate the trailing edge tone noises arising 
from large scale turbulence. Good results in frequency 
and SPL are obtained and compare favorably to the ex- 
perimental ones and the results by other numerical ap- 
proaches. 

In the DNS type approach, we attempt to simulate only 
the large scale turbulence, which is believed to be the 
tone noise generating source. The smaller scale turbu- 
lences (higher frequencies and wave numbers) are practi- 
cally ignored without a subgrid scale model, since we are 
interested mostly in the aeroacoustic data. The CE/SE 
is a scheme of conservation laws. When the divergence 
theorem (2) is applied to find the cell average, it also 
plays a role as a filter ( such as the one in LES approach). 
High frequency and high wave number disturbances are 
filtered out and aliasing errors are unlikeky to occur. 

The novel NRBC based on flux balance is simple and 
genuinely multi-dimensional. It allows the use of a small 
computational domain in the noise computation. 
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